function x = sollu(A,b)

	L = [];
	U = []; 
	x = [];
	[L,U] = lu(A)
	y = soltrinf(L,b)
	x = soltrsup(U,y)
	x'
	 
endfunction
